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Abstract 

In  this  paper,  we  present  a  frequency  jump  detector  for  space  clocks  based  on  two  steps:  a 
time-varying  configurable  deterministic  trend,  and  a  filter  on  the  noises  to  improve  the 
sensitivity. 

We  prove  and  discuss  the  validity  of  the  method  by  analyzing  numerical  simulations 
representing  typical  GNSS  clocks. 


1.  INTRODUCTION 

In  a  Global  Navigation  Satellite  System  (GNSS),  the  user  position  is  estimated  from  the  time  of  flight  of 
the  signals  generated  by  the  system  satellites.  Any  anomaly  of  the  atomic  clock  behavior  on  board  the 
satellites  causes  errors  in  the  localization,  and  clock  anomalies  may  arise  due  to  different  reasons.  It  is, 
therefore,  fundamental  to  detect  the  clock  anomalies  as  soon  as  possible  and  to  inform  the  users.  At  the 
same  time,  it  is  important  to  limit  the  number  of  false  alarms,  which  are  mainly  caused  by  the  natural 
random  fluctuations  of  the  atomic  clock  noise  or  by  noise  of  the  measurement  system. 

In  this  paper,  we  present  a  frequency  jump  detector  for  space  clocks.  Frequency  jumps  represent  in  fact 
one  of  the  most  common  anomalies  in  atomic  clocks  on  board  satellites.  In  addition,  clocks  on  board  are 
measured  by  a  complex  measurement  system  that  may  suffer  from  missing  data,  outliers,  additive  noises, 
and  fluctuations.  An  efficient  detector  needs  to  deal  with  all  these  issues  aiming  to  distinguish  between 
genuine  frequency  jumps  and  false  alarms. 

The  proposed  detector  is  based  on  two  steps.  First,  a  time-varying  configurable  deterministic  trend  is 
estimated  and  removed  from  the  clock  frequency  data.  Second,  a  filter  on  the  noises  is  applied,  aiming  to 
improve  the  sensitivity  of  the  detector  at  the  level  of  the  frequency  jumps  under  test.  Then  the  filtered 
frequency  data  are  compared  to  a  given  threshold  and,  if  the  difference  is  larger,  a  warning  is  raised. 

We  prove  and  discuss  the  validity  of  the  method  by  analyzing  numerical  simulations  and  experimental 
data  from  GNSS  clocks. 
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2.  MODELLING  CLOCK  ANOMALIES 

We  model  clock  noise  in  discrete  time  as 

y[n\  = /u[n]+ %[n] 

where  y[n]  is  the  average  frequency  deviation,  ju[n\  is  a  nonstationary  deterministic  drift,  and  q[/z]  is  a 
white  frequency  noise  (WFN).  We  now  describe  each  of  these  terms  in  detail.  The  average  frequency 
deviation  is  defined  as 

_r  n  1  f"r» 
y[n]=—  y(t')dt' 
r0  J(;i-l)z-0 

where  y(t)  is  the  normalized  frequency  deviation.  The  white  frequency  noise  c\_n\  is  a  Gaussian 
random  process  with  zero  mean  and  autocorrelation  function  given  by 

r£  ta  ^2  ]  =  kk  ]] = (y2s[p -n2] 

where  E  is  the  expected  value  and  S\n\  is  the  discrete-time  delta  function,  defined  as  A[h]  =  1  when 
n  =  0  and  £>[/?]  =  0  when  nT  0.  Finally,  when  the  clock  is  behaving  according  to  the  specifications, 
that  is,  in  nominal  mode,  we  model  the  deterministic  drift  ju[n\  as 

ju\tt\  =  aQ  [«]  +  ax  \nfi  +  A  sin(2  7tf0n  +  (p)  ( 1 ) 

where  fQ  is  a  known  value.  The  functions  c/0[n], «,[/?]  are  slowly  varying,  and  we  assume  that  they 
are  approximately  constant  on  a  window  of  Nw  samples,  that  is 

a0[n']~ao[n] 


when  n'=n  —  Nw  + 1 

When  the  clock  experiences  a  sudden  frequency  jump  with  magnitude  A aQ  at  time  na ,  we  replace  the 
parameter  aQ[n\  by 

ao '  \n  ] =  ao  [n  ]  +  [n  -  na  ] 

where  u[n\  is  the  discrete -time  step  function,  defined  as  u[n ]  =  0  for  n  <  0 ,  and  u[n\  =  1  for  n  >  0 . 

When  a  slow  frequency  jump  with  amplitude  A a0  and  duration  An  occurs  at  time  na ,  we  model  the 
parameter  aQ[n\  by 
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aQ'  [«]  =  a0[n\+  (n  -  na )n[n  -  na  ] -  (n  -  na  -  A n)u[n  - 1 

An  An 


3.  THE  ANOMALY  DETECTOR 

We  give  the  formal  description  of  the  developed  anomaly  detector  for  a  real-time  implementation.  If  n  is 
the  current  time  instant,  then  the  detector  operates  on  the  Nw  =  Nwl  +  Nw2  samples  of  the  average 
frequency  deviation  y\n\  in  the  time  interval 


n'=n  —  Nw  +  l,...,n . 

First,  the  estimates  dQ,d{,A,(p  of  the  parameters  defined  in  Eq.  (1)  are  obtained  from  the  Nwl  samples  of 
y[n']  in  the  time  interval 

n'=n-Nw  +  l,...,n-Nw2. 

These  estimates  are  obtained  by  solving  the  least-squares  minimization  problem  given  by 

n~Nw  2 

fl0 ,  a,  ,  5,  C  =  arg  min  ^  \a{)  +  axn'+B  s in  (2zr/n«')  +  C  cos(2 7rf0n')  -  y [n'J2 

a0,a{,B,C  n'=n—Nw  +] 

where  the  identity 

A  sin(2  7rf0n  +  cp)  =  B  sin(2  7uf0n)+ C  cos(2^/0n) 

is  used  so  that  the  minimization  problem  is  linear  in  the  parameters  a0,av  B,C .  We  point  out  that  the 
parameters  A,  cp  can  be  recovered  by  noting  that 

Acos(^)  =  B 
A  sin(<^)  =  C 

It  is  now  straightforward  to  write  the  estimated  drift  as 

//[/?']  =  aQ  +  dpt'+B  s\n{27ifun)+  C  cos(2  7rfQn') 


for  n'=n-  Nw  + 1, ....  n  -  Nw2 . 

Then,  we  extrapolate  the  estimated  drift  in  the  time  interval  n'=  n  —  Nw2+l,...,n ,  and  we  subtract  it 
from  the  average  frequency  deviation 

y0[n']=  y[n']- fi[n']. 
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Finally,  we  smooth  the  obtained  detrended  average  frequency  by  a  moving  average  evaluated  over 
samples 


n'+N„ 3  -1 

Z  voM 

k=n' 


N 


w3 


for  n'=n  —  Nw2  + 1, . . . ,  n  —  Nw 3  + 1 .  We  detect  an  anomaly  at  time  ha  if 


>5 


th 


where  ,s'th  is  a  configurable  threshold.  We  note  that  if  the  anomaly  occurs  at  time  na ,  then  the  detection 
lag  is  given  by  /  =  ha  -  na .  In  a  forthcoming  publication,  we  will  discuss  the  choice  of  the 
threshold  s±  and  of  the  window  lengths  Nwl,Nw2,Nw3. 


4.  RESULTS 

In  Fig.  1,  we  show  N  =  1000  samples  of  the  average  frequency  deviation  of  a  clock  experiencing  three 
anomalies.  The  parameters  of  the  deterministic  trend,  Eq.  (1),  are  a0  =  2,  a,  =8x10  3 ,  A  =  0.5, 
tp  =  0,  and  fQ  =1/250.  The  (sudden)  frequency  jumps  are  located  at  na  =279,  569,  849,  and  the 
corresponding  magnitudes  are  A a0  =  2. 1,  2.5,  3  .  Finally,  the  variance  of  the  white  frequency  noise 
J;[n]  is  <j2  =  1 . 


Fig.  1 .  Simulated  clock  noise  (blue  line)  and  detected  anomalies  (red  lines). 
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The  detector  successfully  identifies  the  three  anomalies,  and  the  detection  times  are  na  =293,  581, 
850  .  The  corresponding  detection  lags  are  /  =  14 ,  12 ,  1 .  Generally,  the  detection  lags  decrease  as  the 
magnitude  of  the  frequency  jumps  increase. 

The  detector  assumes  that  large  outliers  have  already  been  removed  by  a  preprocessing  algorithm,  such  as 
a  sigma  filter,  but  it  is  nonetheless  robust  to  those  minor  outliers  which  are  not  easily  removed  by  these 
procedures.  The  simulated  data  contain  four  such  outliers,  which  are  clearly  visible  in  Fig.  1.  The 
obtained  results  prove  that  our  detector  operates  correctly  in  presence  of  these  outliers. 

In  addition,  the  simulation  results  in  Fig.  1  show  three  sets  of  missing  data.  The  developed  detector  can 
be  straightforwardly  extended  to  the  case  of  missing  data,  which  is  common  with  space  clocks.  We  will 
give  this  extension  of  our  method  in  a  future  publication. 


5.  CONCLUSION 

We  have  presented  a  method  for  the  detection  of  anomalies  in  atomic  clocks.  The  method  is  designed  for 
space  clocks,  where  frequency  jumps  are  common  anomalies,  and  where  the  clock  noise  can  be 
effectively  approximated  by  a  white  frequency  noise  plus  a  time-varying  deterministic  trend  made  by 
linear  and  sinusoidal  components.  The  method  is  robust  and  can  operate  with  missing  data  and  outliers. 
We  prove  the  validity  of  the  method  with  simulated  data. 
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